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1. Inleiding 


De aanleiding tot dit onderzoek was onze gebrekkige kennis van 
de wisselwerking van de drie soorten magnetohydrodynamische (MHD-) 
golven, Een beter begrip hiervan is van belang in de theoretische as= 
trofysica. Het blijkt (zie 1e interne verslag hierover) dat het hier 
gaat om een niet-lineair verschijnsel. De hierbij optredende proble- 
men maken het nodig, het algemene niet-lineaire gedrag van de MHD- 
vergelijkingen te bestuderen, zonder zich tot golven te beperken, 

In dit onderzoek wordt geprobeerd aan te geven hoe men deze ver- 
gelijkingen zou moeten oplossen, Het verslag is als volgt samenge- 
steld, Ten eerste een overzicht van het lineaire gedrag van de MHD- 
vergelijkingen (8 2), Ten tweede zullen we enkele methoden age 
niet-lineaire vergelijkingen op te lossen beschouwen en gaan we na,” 
wat er elders al aan gedaan is (8 3), Ten derde een overzicht van het 
eigen werk (8 kh) en ten vierde een beschouwing van de mogelijkheden 


tot voortzetting van dit onderzoek (8 5). 


Het onderzoek naar koppeling van golven dat elders wordt gedaan 
is meestal gericht op het verklaren van type II en III radiostoten of 
hogere harmonischen in een plasma. In die gevallen gaat het echter om 
frekwenties die hoog zijn teo.v. die van MHD=-golven, Omdat de MHD= 
vergelijkingen alleen toegepast mogen worden als in het medium geen 
plaatselijke ladingsscheiding optreedt (plasmagolven), zullen de re- 
sultaten slechts gelden voor Laagfrekwente veranderingen in een plas- 
ma, Vooral in het interstellaire gës en in de korona zijn de kondities 
hiervoor gunstig; in mindere mate ook voor de aardse magnetosfeer (zie 
8 2). Bijv. het niet-lineaire gedrag van MHD-golven, opgewekt door 
granulatie of spikules, die overgaan in schokgolven zou een mogelijke 


toepassing kunnen zijn, 


2, MHD=vergelijkingen 


nn en en en en en 


Gegeven een plasma met druk p‚ äichtheid p‚ magnetisch veld Ï, 


* 


et & 7 jn 


2. 


elektrisch veld E, snelheid v en stroomdichtheid En Alle grootheden 
zijn in Gauss-eenheden bepaald op een vast punt r van een inertiaal- 
stelsel in het plasma ten tijde t. Dan gelden de volgende vergelij= 
kingen (zie van Kampen en Felderhof, hst. II): 
Te, Kontinuïteitsvergelijking 
B e 
dt 5, Y . (pv) 
2e, Derde wet van Newton 

dv 3 

Pr == E 

waarin F de kracht die per volumeëenheid op een plasma=-elementje met 
dichtheid p en met snelheid v werkt. De afgeleide - is substantieel, 
dew,z. die moet men nemen in een stelsel dat met het elementje meebe= 


weegt, Transformeren naar ons inertiaalstelsel geef t 


dy n av > 
at pt W.W v 


(zie Landau en Lifschitz, hst, I), De kracht F kan men splitsen in een 
inwendig deel ten gevolge van p en een uitwendig deel ten gevolge van 
de Lorentzkracht (we verwaarlozen zwaartekracht, viskositcit, en elek- 
trostatische krachten ten gevolge van ladingsscheiding), De vergelij- 


king wordt dan 


SJ 


=p BE Pp GV) v + Tp ek 3 AE 
fe) c 
je. De wetten van Maxwell 
hi ee ca Ee en 
AH ne Bkk: V.E 0) 
V zj we ON emme iz 
AE ok V.E ber) 
Daar we geen ladingsscheiding beschouwen is de ladingsdichtheid ĳj = O0, 


Dit betekent dat we ons tot laagfrekwente veranderingen beperken, zodat 


we mogen stellen: 


Hierin is og de geleidingstensor, We nemen aan dat het plasma isotroop 
is; dan is co een konstante, Deze wet geldt als we met het plasma-ele- 


mentje meebewegen; transformeren. we naar het inertiaalstelsel, dan 


vinden we voor v << c 


vaAE) 


[a 


j = 0 (E + 
Als co = 0, moet gelden dat 


E + v A Ë == 0 


en EN 


Zodoende blijft j eindig; we verwaarlozen uitwendige E-velden en Hall- 
effekt. 
De, De toestandsvergelijking. 

We nemen aan, dat p eenduidig door p wordt bepaald. De toestands- 


vergelijking luidt dan 
p = fp) 


Voor p kunnen we ook de differentiaalvorm schrijven 


- Se 4 
dp ao Ce 
Ee, 8 dp dp dp 
waarmee de tijdsafgeleide wordt = 


Nu transformeren we naar het inertiaalstelsel en substitueren de kon= 


tinuïteitsvergelijking: 


9 > _ ef n 
Er +GU)Pp = B ì Ens + (r.V)o } 


C è 
EE RAENRA: 
Als we te maken hebben met adiabatische veranderingen, geldt 
Pp = Ep! 


waarin v de constante van Poisson en K een konstante die van de ge= 
middelde temperatuur en het gemiddelde molekuulgewicht afhangt, De ge- 


luidsnelheid s wordt in dit geval 


ea Ad 
3 


Deze snelheid hoeft in het algemeen niet konstant te zijn! 





Uitgaande van deze vijf vergelijkingen kunnen we E elimineren met 
de wet van Ohm en j met de eerste Maxwellvergelijking, Dan vinden we 


het vereenvoudigde stelsel 


av en ne ek => 1 = => 
mf SP (v.V) v + Vp + In BA Fn E) (Ja) 


les 


5 5e = Y.Cpv) (Ab) 
“ot == Velp te p Vv 1e 
en = VaA(aAE) (1d) 


Samen met de vergelijking voor s is dit een volledig stelsel. In het 
vervolg zullen we (1) "de MEHD-vergelijking'' noemen. (Een voortreffe- 
lijke diskussie van het feit dat (1) behoudswetten zijn vindt men bij 


Van Kampen en Felderhof, hst. III). 


Geldigheid 


Voor de toepassingen is het belangrijk te weten, in welke om- 
standigheden de gebruikte benaderingen geldig zijn. Beschouw een la= 
dingsdichtheid A over een gebied met karakteristieke lengte L, Dan 


volzt uit 


le wet van Maxwell: d As in J 

L le) 
he wet van Maxwell: E mz UL 
wet van Ohm: jj z eN L 


Hiermee vinden we voor de orden van grootte van de Lorentzkracht resp. 
de elektrostatische kracht: 
BN Pe 
bitte se re nr 
C | 2 
C 
De verhouding van deze krachten is dus 


ie e 
L.k, ( G L ) 


Siel c 


Bij een hoog geleidingsvermogen en/of karakteristicke lengte is de Lo= 
rentzkracht dominant. In speciale gevallen is de elektrostatische 
kracht nog wel van belang, bijv. als het medium zich zó aanpast, dat 
3 EN Ë . Het ligt voor de hand, dat deze situatie na verloop van tijd 
ontstaat omdat het plasma moeilijk loodrecht op de veldlijnen kan be- 
wegen, In dit geval verdwijnt de Lorentzkracht en we hebben wegens de 


le wet van Maxwell 
E … Vai 
Velden die hieraan voldoen noemt men krachtvrij (force-free fields). 
Voor zeer grote L volgt uit de verhouding der krachten, dat we 
Jg oo mogen nemen, Omdat Ea i}L volgt bovendien dat fl vO . 


Tenslotte een relatie voor de toegelaten frekwenties, Uit de 1e 


en be Maxwellvergelijking volgt de kontinuïteitsvergelijking 


9 


sE 
| 
t 
hs 
il 


Voor een karakteristieke frekwentie f volgt hieruit 


EN in We hadden al dea aft 


zodat we vinden f ‚40 . De benadering is dus goed voor f << 0 , en 
we zagen al dat oa als groot beschouwd mag worden als L groot is, Men 


vindt bijv, 








1 EE s° 
L (em) | o (sec ) ( an ) 

H II -gebied | 101 | gole | 40% 
H I= gebied | 1017 | 10° | 10°! 

$ ï 
buitenkorona | 107 101% | 10°? 
binnenkorona i 10° | 101% 10°2 
P 7 8 9 
ionosfeer 10 10 : 10 


(zie Van Kampen en Felderhof, pag. 6). Toepassing van de MHD-vergelij- 


kingen in korona en interstellair gas is dus zeker toegestaan. 


Lineaire vorm 


Om de behandeling van de niet-lineaire MHD-vergelijkingen in $ 4 
te begrijpen moet men precies weten hoe het lineaire geval wordt aan- 
gepakt, Daarom gaan we hier tamelijk uitgebreid op in. 

Elke variabele kan worden ontwikkeld naar een dimensieloze para- 


meter « in de reeks 


VN + ET 


€ Foe v 
1 


V = Et 2 


die in de vergelijkingen kan worden gesubstitueerd, Daar voor elke e 
het stelsel volledig moet blijven, mogen we gelijke orden van … bijeen 


nemen. Zo krijgen we 





ais Re Ee 
| = Lo = V. (pv) 

nulde : sl Ee 

orde S ÔP, oe eN 2 ie 
SA BE Po Te 
ERIN 
ae A (vo A H) 


Als oplossing nemen we Yo ==) 5 Po’ Po en B homogeen en stationair. 
Hier zullen we in het vervolg steeds van uitgaan. Als «& = 0, is 
hiermee het stelsel (1) opgelost, Als O <e << 1 nemen we ook de 


eerste orde er bij en vinden 


r 


LD 


= e= 7 —_ è 
| Po Bt Vo, + Eo A WAE) KE 
2 U (25) 
Se =P Dg Telg 
eerste £ ap 
orde 4 1 2 2 
sne hi 2 
| 3t So Po v.v, (ac) 
| al, 
| TE Oe V A Cv, & H) (ea) 


Dit stelsel noemt men doorgaans (vaak met inbegrip van de zwaar- 
tekracht in 2a) "gelineariseerde MED-vergelijkingen!!, Deze uitdrukking 
is misleidend, want behalve dat linearisatie is toegepast, hebben we 
ook al een aanname over het nulde-orde stelsel gedaan! 


In eerste=orde benadering zien we direkt, dat 2c en 2b identiek 


zijn indien Py = 5 Con ‚ De vergelijkingen zijn dan te reduceren tot 
v 
EA nd = eg Y Bg + TE EN AT A E) (3a) 
òp, stak 
Ge EH Vev, (bo) 
kn Talk aj 5 
Ò atd: nad (sc) 


Nu differentiëren we (3a) naar t en substitueren (3b) en (5e). Dan 


… E … . 3 ad =s 
krijgen we een vergelijking voor v, alleen: 





2 Â 
ò v 2 oe 4 eN > ï > => \ 
5 k is 4 T or 
D zt8 B Verd rn EA (VATA 4 EH) 0 
at 
Noemen we E Clr pn = C , dan vinden we 
0, 22, => ê > à 
De en So Vv ev) + Ge Â L 7 A 7 Â (v, Á E) ES = OQ 
Substitueren we voor Vv, een vlakke golf: Ee ge mn 
dan vinden we (zie appendix kh}: 
IN + (5 40 OETDR + CD {ELD - EE ERE) 
-1 o a an a” Ì a“ 1 a“ 1 * 1 a) 
= OQ 


Er zijn nu drie gevallen (“'modes"', polarisaties) mogelijk: 


fe) > > — > 
1e ke 18 CE, ‚ k UN, 


Dan hebben we een longitudinale magneto-akoestische golf met dis- 


persievergelijking ge -# (as + s°) = O0 
De fasesnelheid is u= 2 = Wi C Ë + 5 Ee 
k a o 
en de groepssnelheid ís Se = U, Deze snelheden zijn gelijk, dus 


er is geen dispersie, Daarom spreken we in het vervolg van !de!' 


ij 


tes 


RAS Lite Deze redenering geldt ook voor de overige modes, maar 
niet meer als er dissipatie optreedt (zie bijv. Jackson, pag. 
BR Bada 

2, EMG, ENT 
Dan hebben we een longitudinale geluidsgolf: met Aieperelervengde 
lijking de Ks = 0 
en snelheid B . 

S. EAD, ELF, 
Dan hebben we een Waar Tersdle nan a met dispersievergelij- 
king Res ED & = O0 


en snelheid C, 5 


5. Niet-lineair gedrag 


De aannamen uit & 2 zijn zo beperkend, dat we als oplossing li-= 
neaire golven konden vinden. Aan de golfvergelijking zien we al, dat 
zelfs een kleine niet-lineaire term grote algebraïsche problemen 
geeft, Om een indruk te krijgen van de resultaten die we in het alge 
mene geval verwachten, lossen we (1) op in het eendimensionale ge val 
voor È = 0. (Riemann, 1860; zie Landau en Lifschitz Ô 94). Deze be- 
handeling berust op de methode der karakteristieken (zic appendix 1), 


De vergelijkingen worden 


3v ELSE … 

ot * Y 5x t px 8 

2e , pv) 0 

ot Ox u 
MN: ä n 3(ev) = : Se Jpv = ER EN nd cl 2v 
u weten we, dat Dx v Bx + pp DX ii do 3x tg dv do Bx 


Omdat v en p beide funkties zijn van x en t,‚ mag v als funktie van p 
worden beschouwd, zodat 


der) _ „dep, , dra _ lov) ap 
Ax dp 3x dp 3x dx dx 


Een analoge redenering geldt voor Pp, zodat 


2e , dlev) A _ © 
dt dp ax 


WW vad hd 


av 
bt D EE ox 8 


Omdat de differentiaalvormen gesloten moeten zijn, geldt 


en 


PE _ 8E SE Om tE 
dt 3p dt 8 at Ar ot ee 


Ĳ 





8. 


Daar v een funktie van p is, zijn beide uitdrukkingen gelijk en vinden 


we 


BE 5 LB 
Vv + @ de Vv + a Ev 
Stellen we dp = Er do = a do , dan volgt na integratie 
vz + E 5 dp + JE 


De konstanten J2 zijn de zgn, Riemann-invarianten. Deze zijn konstant 
op oppervlakken Xx = + st , d.w.z. op karakteristieken, Substitueren 
ve de gevonden v, dan blijkt 


@ 


ÒX 5 
der: Sn = v+s(v) 


ali 


hetgeen geïntegreerd kan worden: 


XxX = t Ì v + s(v) } + f(v) 


waarin f een willekeurige funktie van v. Door inversie van f vinden we 
elf 
voe feet (v + s(v)) 
De voortplantingssnelheid s van de verstoring beïnvloedt de gassnel- 
heid v; als f strikt monotoon is, is v voor beweging in de voortplan- 
tingsrichting groter dan in tegengestelde richting, Dus een golftop 


gaat sneller dan een golfdal, zodat de karakteristieken elkaar gaan 


snijden: 


Dat betekent, dat v niet meer eendui- 
dig bepaald is binnen het aangegeven 
gebied, Er ontstaat dan een schok- 
golf, 
Uit het bovenstaande kunnen we 
de volgende konklusies trekken: 
i) Door niet-lineaire effekten kun- 
nen instabiliteiten ontstaan, We 


moeten er dus rekening mee hou- 


ge el den dat alle oplossingen voor 


eindige tijd geldig zijn. 
ij) Verandering van een lineaire golf is mogelijk, bijv. doordat s 
varieert ten gevolge van een dichtheidsgradiënt of een scheidings- 
vlak tussen twee media (zie literatuurlijst). Hiermee zullen we 
ons niet verder bezighouden. 
iid) De moeite die we hebben met het oplossen van de eenvoudige inter 


aktie van een golf met zichzelf doet ons het ergste vrezen voor 


9. 
meerdimensionale gevallen, De geldigheid van de resultaten is in 


elk geval beperkt tot de tijd oe waarop schokvorming optreedt. 


Dit gebeurt als v over een afstandje diskontinu verandert, dew.z, 


(EE) = @ > (& ze Haan 
opd t dv t o= dv ov 
o o 
ofwel 
kt mm DET De 
o Ov Vv 
we zien dat Ln groot is voor Vv <s 3 in dit gebied zullen de 


resultaten dus goed zijn. 


De tweede en hogere orden in « hebben wiskundig dezelfde vorm 
(zie } &)., Dus is de transformatie die aan de MHD-vergelijkingen de 
diverse orde-benaderingen toevoegt, niet-singulier. Verder is het 
eerste-orde probleem oplosbaar, Volgens de stelling over impliciete 
funkties zijn dan ook alle hogere orden oplosbaar, We zoeken dus niet 
voor niets, en we mogen ons tot de tweede orde beperken, omdat hogere 
orden met dezelfde methoden aangepakt kunnen worden. 

Het probleem van iedere niet-lincaire benadering is het volgende, 
Een afwijking is per definitie lineair, als zij klein is ten opzichte 
van de eigenlijke grootheid, Dat wil zeggen, dat de toestand waarvan 
wordt afgeweken steeds herkenbaar blijft. Afwijkingen die hier niet 
aan voldoen noemen we niet-lineair, Dus hierbij herkennen we de oor= 
spronkelijke toestand niet meer! Daarom helpt het ons niet, of deze 
toestand ‘netjes! is (bijv. golfvormig), want we vinden er niets van 
terug, Vandaar de opmerking in & 1, dat we verplicht zijn het algemene 
geval te behandelen in plaats van ons tot golven te beperken, Ook een 
“nette! geometrie biedt geen voordelen, behalve voor de randvoorwaar- 
den en het aantal vergelijkingen. Een werkelijke vereenvoudiging ont= 
staat, als we ons beperken tot minder dimensies. Bij aanwezigheid van 
een magnetisch veld heeft men altijd nog twee dimensies nodig, en het 
weglaten van En en En is twijfelachtig omdat dan de magneto-akoes- 
tische modus verdwijnt, Daarom gaan we uit van het volledige tweede 


orde probleem: 


10, 


í ar eh 
| mn =D, +Vo, er B, A AAE) « 
| zE, A Fai) (ha) 
dh | ee = Y. (pvo + pa) Cb) 
= oe ee v Yp, + 5 °p Ft, + es adr, + ÔR BP, Vv, Che) 
ES = Vaal) +TaAG, AE) | (4d) 


In het volgende deel van deze paragraaf sommen we enkele elders ge- 


volgde benaderingen op. 


Üquivalent source method 


Bij de verklaring van type II en type III radiostoten gaat men 
uit van een volkomen willekeurige verdeling inhomogeniteiten in het 
ges, die men als oplossing van (2) stelt, De fouriergetransformeerde 
bovat de eigenfrekwenties in ruimte en tijd van deze verdeling, en de= 
ze frekwenties wekken het Vo-veld op (equivalent source, “schijnbare 
bron!, Zie Ginzburg en Zheleznyakov (1958, 1959, 1961) en Sturrock 
C.S. (1965), 


Matrixformulering 


Als men de variabelen samenvat met de vektor g = (v‚o‚p,l) dan 
is het stelsel (1) formeel te schrijven als 


fo) => => 
En = ÀÂ (g,v) 


ein! 


Men kan proberen oplossingen te vinden door achtereenvolgende benade- 
ringen van A, door gebruik te maken van symmetrieën (Frisch, 1964), 
door aan te nemen dat A monotoon van E afhangt en dat ee weinig vari= 
eert over de lengte van het gebied dat we beschouwen ©& (Kaplan, 1960; 
Frisch, 1964), Deze methode wordt vooral gebruikt bij de berekening van 
golfgedrag ten gevolge van dichtheidsgradienten (zie literatuur). He- 
laas werd in geen enkel artikel bewezen, dat het berekende gedrag een 


speciaal geval is van een algemeen te verwachten koppeling, 


mn en nm nn mn mn ne a a mn en TE 


Als een golf een variabel medium doorloopt, zijn er twee grensge- 
vallen: 
le, Het medium verandert merkbaar over een afstand, die groot is verge- 
leken met de golflengte; de energie van de golf is vergelijkbaar 


met de potentiaal die het storende medium beschrijft. 


gee 


2e. Het medium verandert willekeurig, en de energie van de golf is 
zroot ten opzichte van de potentiaal. 

In het eerste geval kunnen we de WKB-benadering, ook wel “eikonaal- 

methode!’ genoemd, toepassen; in het tweede geval de Bornbenadering ee 

(zie appendix 5 en 2), Benadering 2e. werd in geen enkel artikel ge- 

bruikt, en Je. alleen bij Fitremann (1966) en Weinberg (1962), 


Hamiltoniaan-methode 


Men kan bewijzen (Van Kampen en Felderhof, pag. 71) dat (1) met 


een variatieprincipe kan worden afgeleid uit de lagrangiaan: 


ke 


ae =2 J.R a 
dl Ko dik Co) - 8n p ee Bk 


_L(v,r) = 
waarin n de dichtheid op t = 0 ter plaatse Xx, en Y de kompressie- ze 
energie per massa-eenheid: 

| 
Elp) = | de áp 


0 


Daar we geen potentiële energie invoeren (bijv. door zwaartekracht of 
ladingsscheiding) is L gelijk aan de hamiltoniaan H. Deze kunnen we 


proberen te schrijven als 


HE == M 
1 L + H + Ho 


waarin H, en H de ongestoorde golven voortbrengen en Hip de interak- 
tie, Dit is een storingsprobleem met standaard oplossingsmethoden (zie 
Merzbacher) , 

Een andere mogelijkheid is, een Fouriertransformatie van H te ma- 


ken. Als H e= H(g) , substitueren we 





ene: de teen ee ER ab 
g(r,t) = 5 | a ER du g (k‚o) e= 7 TL 
J 
(en) 
Men kan = opvatten als kreatie=operator voor vlakke golven, en z als 


annihilatie=-operator. Wegens de eis van kausaliteit moet dan 
A Se 


== t_ voor elk produkt van E en g 


A 
ne 


ECR, vw).8 


% 


E's") 


en wegens energiebehoud moet k, = ES ‚ De werking van dit produkt 
op een toestand met twee vlakke golven wordt schematisch weergegeven 


door 
>= — 


k k er N 
1 1 Op het tijdstip t dat de operatoren 


Ar verdrag 


werken is de golf KE, in CA ‚ Op ana= 


behinds a loge wijze stellen we de werking van 
2 2 


SE) gE, vo) E (E,0) voor door 


K Ei 
Lt PAL wr AD ad 
Ei En 
Aes 
da 
d.w.z, annihilatie van de golf K en kreatie van K, en es let deze 


operatoren brengen we een lineaire golf voort door 
Ax aen 
Ee Bieger ge à 


(zie Abrikosov). Een niet-lineaire golf ontstaat door termen van de 
gedaante 


AAA A A 
de 55 5 Ar ae) 


EEE + 56E + 8 EE +EEBEH oee 
waarbij steeds één golf verdwijnt en er twee ontstaan, of twee ver- 
dwijnen en één ontstaat. De verwachtingswaarde van a} deze processen 
samen geeft aan, of er een netto overgang is van de ene golf naar de 
andere, Deze methode werd alleen toegepast door Krishan en Selim 


(1968) voor koppeling van longitudinale en transversale plasmagolven. 


L. Onderzoek 


Matrixmethode 


On am nn nm an en en en a ve 


als 
DE ET EN 
SE == Als, We 


(zie appendix 5), We kunnen A zodanig splitsen dat 


ke K(g) ES L(g) ar * (a) En 


waarbij K, L en M geen differentiaties meer bevatten. Uit deze regel 


leiden we de golfvergelijking af door differentiëren naar t: 


d eà Á 2 2 2 EE pen - ge 
5 eK 2E ‚rf 2E Mm 2E En + Es Ei LE + nk ) zE + 
at Ox 3y dz 7 
aL aL aL aL © DE 
amen mn M= 
+ ge % Ke sk len + MeT ) by + 
g Ee 
AM _ aM aM AM 9E 3 
Sr KES IE 4 MES ) LE Ei Ab 
+ ( et tlr t 57 * ez ) iki (KL + LK) 5705 * 
AT ie 
+ (LM + ML) SE 4 (MK + KM) ê& 
9y9z ÒZX 


In één dimensie is dit 


=> 2 => 
ÄE Es Ke LE + ge + K 5 28 
Bte ax dt Ox DX 


hetgeen voor konstante K de bekende golfvergelijking geeft. Als K 


langzaam met t en x verandert, wordt dit 


2 2 — 
ne = KE 2E, c 28 
ot Ox ea 
waarin « een storingsoperator. Indien 5 = EN eAeriot ‚ vinden we 
B Rn nde 
„RE Ss ek Ll mee tf En 


Omdat we niet weten of k een inverse heeft, lossen we deze vergelij- 


king op voor # ke 8 


De DE en 2 2 
Janes 8) = kileg ) « Ck B) "08 (k B’ = 
= wé ‚p … 164 A we) lend pn ie el 
L Ak Bo ” Meg, } + eg) D 8e (k Bo) O 
; 2 > 2 2 1% 
LE Es = 4 n be 
waaruit volgt #k 50 fides, 218 (k B) j 


2) 
G CES « .. ‘ 
Termen met ce” zijn verwaarloosd. We zien dus dat K een reëel deel 


heeft dat de gewone golf beschrijft en een komplex deel dat de ampli-= 


tude verandert: 


2 ep 
a = 4 
Im ((x Bo) hes, 
Hieruit volgt dat de golf gedempt wordt voor € > O0 en versterkt voor 
« < 0 „ De karakteristieke lengte Xo (d.w.z, de e-waarde lengte) voor 


beide effekten volgt uit 


ex = el Sd 
© 


Zoals te verwachten, is Xo groter naarmate e kleiner is, Voor een 


tweedimensionaal systeem is 


2 2 n E > 2 
= = K£ ae + ( 2 + rd + ze ) 28 + Ee 2E + 
ot Bo 6) fop;< òy Ox ay 

> 2 

„OK „ ôL aL òg os 

En, KZ 2e le 

dt rid mak L 5y ) Ee + (KL + LK) Bor 


Blijkbaar is ook bij vaste K en L een gekombineerde oplossing mogelijk 
(magneto-akoestische oplossing). Niet-lineaire koppeling krijgen we 
pas, als L van x afhangt en/of K van Y. Voor kleine variaties van K en 


L krijgen we 


25 KLEI LE, (mn) DE, E‚n ZE 
ge Bx oy OXY Ax ôy 


Het gedrag van de eendimensionale oplossingen (K- en L-golf) is zoals 


Th, 


boven beschreven; de KL=golf wordt gedempt voor « + ùÛ > 0 en ver- 
sterkt voor ze + <0 . Hieruit volgt, dat diverse soorten energie= 
overdracht mogelijk zijn, maar dat nooit de KL-golf gedempt wordt als 
beide andere versterkt worden, en nooit versterkt als beide andere ge 


deupt worden, De karakteristieke Lengte x, voor de L-golf volgt uit 


1 
” 2 
xj = 2L en voor de KL-golf hebben we 
‚ è 
X Ce +1) = 2 (K + L) ofwel 
Cx, Ln ze + , En X, DM = 2(KL + LK) 


De verandering in de KL-golf hangt direkt samen met die van de K- en 
L-golven, Uit de vorm van e en 1 (d.w.z. van K en L) moet de energie= 


overdracht worden berekend, Indien bijv. de KL-golf niet verandert, is 


st == OD 4 zodat 
ACK+L) … BERKEL) ACK+L) 
ae aad 
en e == , d.w.z. ns DE ofwel ra 5 Ke 5 
waaruit volgt K+L = 0 of K=-L = O . Het eerste geval is 


een oplossing van de differentiaalvergelijking; in het tweede geval is 


BK BK ak 
ot dd rk 


Voor alle K die hieraan voldoen kunnen de K- en de L-golf onderling 
koppelen zonder de KL-golf te storen, 

Bedenk, dat en s 5 niet exakt behoeft te gelden, omdat Xo en 
Xs schattingen zijn. We weten alleen, dat resonantie optreedt als de 
Lensten Xx, en Xx, gelijk zijn. 

Het bovenstaande is kwalitatief, en betere oplossingen zijn moge- 
lijk. Ten eerste kan men proberen de vergelijkingen op te vatten als 
partiële cifferentiaalvergelijkingen met veranderlijke koëfficiënten, 
Om dit te kunnen doen, moeten we een veronderstelling doen over de op- 


lossingen en 2 = E(r,‚t) schatten, zodat we kunnen schrijven 
AC5,V) mm A(F,t,V) 
De vergelijkingen kunnen dan met Born- of WKB-benaderingen behandeld 


worden, Ten tweede kan de vergelijking met suksessieve benaderingen 


worden aangepakt, Bij voorbeeld, in één dimensie hebben we 


2 2 8 — 
Bee Ee ot ops DX 


Stellen we K even konstant, dan krijgen we voor E een vlakke golf 
B irstde Hiermee berekenen we een benadering K van Kk, en krijgen dan 


1) 


1D 
dn 2» ak ok ze 
EE su dBi Se 
de o a at o Ox 9X 


=> u . * EJ 
Met de oplossing g hiervan berekenen we K, en zo kan in principe een 
algemene oplossing worden gevonden. Tenzij K niet kontinu is, verwach= 


ten we konvergentie, 


Green=-funktie methode 


Ve gaan uit van het stelsel (4), en schrijven dit als 


ov : 

SG 2 KS: > _ ee > En 
ist ® Vp, kens EH, A (V A H_) = Á (5a) 
rk Fw Vo = B (5b) 
Op x 

2 Br ve 
5E + B Po Ve Vo = D (5e) 
5 M En T 
EE VA Cr A B) N (5d) 
Ov 
waarin -Â = p dm te EH A (V AE) 
i 1 at o 15 1 in 71 1 
„B bn Ve Coov) 
D = vt ss Vv, +2 Vv 
To Sp Pr 4 Po Te 
nf 5 TE, ad —> 
Vv Á Gr, A H) 


Het stelscl (5) heeft de vorm van (2), maar nu inhomogeen. De inhomo- 
gene termen volgen uit het eerste-orde probleem, We zien dat elke ho- 
gere orde in deze inhomogene vorm geschreven kan worden, waarbij de 
inhomogene termen uit de lagere orden volgen, Daarom beperken we ons 
tot het tweede-orde probleem, Differentiatie van (5a) naar t en sub= 


stitutie van (5e).en (5d) levert 


0°v, 2 2 > > => > => ij > 
es ST (Hv) + D A{TAT AG Er) et 
Bee o 2 El è a J 
waarin  & EN Baie de dE 
Pz a Ps at 


Formeel. kunnen we deze vergelijking schrijven als 


L Vo = 1 


waarin L een operator is, Deze vergelijking lossen we op met de Green- 


funktie, We definiëren eerst de vektor 


x = (F,t) 


Ad 
” 


d id . . Ke 2 . . 3 e - . 
Laat e, een basisvektor uit de r-ruimte zijn, en I de identiteits- 
nd bd > 
operator, Stel ko Post) waarin ren Et gegeven vaste plaats en 


tijd zijn, Dan definiëren we de Green-funktie GGL) door 


ig, 


ï, L G GLX) e, fj = ô (x — x) I e, 


Formeel betekent dit, dat we de oplossing Vo ontwikkelen naar de 


eigenfunkties van het homogene probleem. Deze voldoen nl, aan de 


orthonormalisatie 


rl AËE de sh 
(2) 
Met behulp van de Green-funktie lossen we de inhomogene vergelijking 
op door 
v(%) = | ER J Pe hd 4 randtormen 
è Ĳ o o o 


waarin G een lineaire kombinatie van de gd en YP het gebied waar Ë ge= 
definieerd is, We nemen aan dat T op de rand van Y nul is, zodat de 
randterzen verdwijnen, 

Bedenken we dat G een operator is, dan zien we direkt dat Vo een 


oplossing is: 


> 5 ĳ { sd > } h > je 
Lv, = ) E { GO ie) TG) d Xo > 
y 
E fe sad \ >p bs „ 
= | { LG Gele) ) Tx) d Xs 
ba 
=| OA en Xx) Pe) dx = kN, 
Y 


Om de resultaten te vergelijken met de dispersierelaties uit het li- 


neaire geval, passen we Fouriertransformatie toe: 


B Rn _ d „> ik 1 
5 (x » Xx) ne (2m) h f oik. Cx Es) dT 


al - 


CS er ex) lj 


pr bor Ee 


d k 





akk) = 
© 


HER GC 


waarin E =S Lo) en & het gebied waar golfvektor Î en frekwentie 
gedefinieerd zijn. Substitutie van Ge, geeft (zie appendix k) 
è 


=(ge, Jo 


f — 16 2 E => IK f Wi =d Bl ' Ee EEN 
+ {lee Ake - CAF AL AGS) AED) = (eme, 
de kunmen nu de bekende drie gevallen onderscheiden: 


1, ge, zkt Â en Ë VN ‚ zodat 
_Î 


aal} 
3 (ar = Ce) (A1fvéngolf) 


ge. = ë, (2m) 
2 ge, NR en eN IÌ , zodat 
_Î 


Ee, = B 4d lie fa) (geluidsgolf) 
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En ge, IA en e, | Ì , zodat 
ge, = ë Er Pe _ a B + oe uk (magneto=akoestische golf) 


r 


We hebben hier drie eigenwaardenvergelijkingen voor 8, met verschil- 
lende eigenwaarden, Daar g driedimensionaal is, kunnen we dan de dia- 
gonaalgedaante veronderstellen. Voor de drie oriëntatiemogelijkheden 


van Vi krijgen we dan drie verschillende Green-funkties 


e, GE) _ Î (a) PP (For) Ho lteto) ta 


ö 
Zoals bij deze problemen gebruikelijk, is de integrand singulier, met 
enkelvoudige polen in w = + Ke ‚ Dan vinden we door omschrijven van 


de integratie over 7 
en Die) èn TL 
f : he € : En , 
da RIAS Pr) a = | (Ê af | aÛ | Ei ò do Mir ro leos t) 


öt o o o 


waarin &' het gebied waar 4 gedefinieerd is. We vinden 





r En Pee 8 ed 
Es al 1 LE Es) zE Pol) a 
EER 
Co Ee ai ed 
= EE 1 rel „reel op … 
had d > KS 
ilr-ro | 5 olr-r, 
a > > 4 
2 E 7 Er | en de df . 
rr, | olr-r.l 5 


Hiermee vinden we drie typen Green-funkties (zie appendix 6): 
oo oo 


” i ê r 1 
) za a | a -e. 


(em)? [r-r, olr-r, En Ee 
=00 =D 


gi „ 
iki r-ro| + dolt-t) af 

















[F-r, | 
oo dn KA Ek 

-iIl 9 Û 1 Luk B he. 

vetes OO a 
(en)° [r- -rol ol r- CUE ne 
En [r- =| 
NE EEN Oe 
7 2n{r-r | £Ca 8 ä 
| o 65 
Ir-r | 

E i 8 

pad 6 ( Ca + in ad np, e 


hr Ca |F-r.| 


Deze funktie bestaat voor alle t < to en beschrijft dus een inkomende 


golf (geretardeerde oplossing): 


18. 





[r-ro | We zien duidelijk de overeen- 
ee Pi B iens komst met de lichtkegel, en 
G, is hier de propagator van 
de Alfvénrolf,. 
4 
P. “Fer, | 5 0 (beeb 
Op analoge wijze vinden we: 
-|F-E 
Brin % I | 
G, GEE) nen ( hand en + te t, ) 
in C_ |rer a 
a o: 


Deze funktie bestaat voor t > an en beschrijft een uitgaande golf 


(geavanceerde oplossing): 


Voor 








he 
= af J du cos py (Cm 
Ee a a 


Deze uitdrukking bestaat alleen bij r= o en t = Ge ‚ en heeft 


4 beschrijft een staande 


o sl 


daar een pool van de eerste orde. Deze G 
golf, 
Op dezelfde manier vinden we de Green-funkties bij de andere twee 


polarisatierichtingen van Vd ; 





5 jn I ee [r-r,| 
6 Gel a) = En 6 C+ rn + tet } 
ln so |r-r,| o 
G Glx ) = PED Ten 6 C+ a + Et - tj} 
5 8 lim (aw) | F-rol ear ai 
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oi E 7 eer, | 
G, Gel) = sn | do eos W Mae” +t=t,) 
2n Ss, r-r,! -ò o 
o > IT ie Adan | 
Gs EI me 2 pen | dj COS u ener + te to) 
EO en (C2ss OFTE | ® (oan Ey 
Oo‘ =00 a o 


Met behulp van de gevonden Groen-funkties kunnen we, zoals ge= 


KN . ad . * . 
zegd, oplossingen v, konstrueren bij alle integreerbare inhomogene 
+ 


_d 


2 
en : en =Û ie & 
termen f ‚ De asymptotische vorm van de CG” is steeds yr „We inte- 
yreren dus de stralingsbijdrage van elk elementje van het medium, 


as 
waarbij 1 als ‘stralingspotentiaal!" wordt opgevat. 


De, Honklusie 


Er werd onderzocht hoe de drie MHD-golftypen elkaar beïnvloeden, 
Dit leidde tot een beschouwing over het niet-lineaire gedrag van de 
uiiDeversgelijkingen, De literatuur hierover is beperkt; meestal houdt 
men zich bezis met plasmagolven in plaats van MHD-golven, Voor het 
berekenen van het niet-lineaire gedrag werden twee algemene methoden 
beschouwd, namelijk oplossen met de matrixformulering en met de 
Green=funktie, Het eerste geval is vergelijkbaar met de aanpak van 
Frisch (1964), het tweede geval met die van Sturrock ces. (1965) en 
van Weitzner (1961), Bij enkele analoge problemen op andere gebieden 
ziet men interessante overeenkomsten (zie literatuur), 

Thans wordt gewerkt aan een uitbreiding van de WKB-methode, in 
verband met het onderzoek van Elly Dekker, Het onderzoek kan op vele 
andere manieren worden voortgezet, Ten eerste het toepassen van de 
algemene resultaten uit 3 & op praktische gevallen, en ten tweede het 
uitbreiden van ê bk tot media waarin zwaartekracht, viskositeit en 
ohmse dissipatie optreden, Dit laatste zal een zware opgove zijn; in 
vrijwel alle gevallen in de literatuur werd iedere vorm van dissipa- 
tie uitgesloten, waardoor toepasbaarheid van de resultaten beperkt 


blijft tot processen op grote schaal (zie 3 2). 


pi 


20. 


APPENDIX 





1. Karakte eristieken 


en en mn mn on Sn en en me vn 


Zen lineaire partiële differentiaalvergelijking voor een funktie 


ò 
u(x) schrijven we als 


n 

5 „u _ 
Jen OE Bx. 5 (1) 
i=1 8 a 


Jr E í > ® . . 
waarin n de dimensie van de x-ruimte en ay en a gegeven kontinu dif 


. . ad . Ea . EN 
ferentieerbare funkties van x . Zoals gewoonlijk gaan we hij het op= 
‘ 


lossen uit van het homogene geval & = 0 . Door het systeem vergelij- 
kingen 

dx, pe 

er. aj (x) 


waarin S een parameter is, wordt een bundel krommen Bj = Hiel 
voortgebracht, Deze noemen we de karakteristieken van (1), Als Xi van 
meer parameters afhangt krijgen we meer soorten karakteristieken, 
Uiteraard heeft deze parametrisering van de koëfficiënten uit (1) 
iets te maken met de richting waarin de oplo:sinz u(X) zich in de 
Xeruimte beweegt, zoals een gewone differentiaalvergelijkine door een 
richtingsveld gedefinieerd wordt. Men kan aantonen (Courant en 
Hilbert II, hst, 2) dat ieder Foplossingsvlak" u(x) door de karakte 
ristieken wordt bepaald; cen karakteristiek geeft dan aan, hoe de op- 
lossing zou lopen als die alleen afhing van de parameter die bij deze 
sarakteristiek hoort, 

Voor het bepalen van een oplossing is het belangrijk te woten 
welke karakteristieken aan de beginvoorwaarden voldoen, Deze karakte= 


ristieken legen de oplossing vast. Bij voorbeeld, als (1) hyperbo- 


lisch is, kan men transformeren tot 


Be 


EE ma PCR) CE) + ar 


waarin P, Q, R bekende funkties van de karakteristieken-parameters À 
en me. Deze vergelijking heeft golfoplossingen voor Y . De karakteris- 
tieken zijn golffronten, en Y moet worden aangepast bij richtings en 
vorm van deze fronten bij gegeven randvoorwaarden (zie Hargenau en 


Hurphy pag. 683). 


2. Bornbonadering 


Gegeven een inhomogene golfvergelijking 


0 
CT 4 3 Y == af 


Fi 


en 


waarin « een funktie van r . Deze heeft als asymptotische oplossing 
een vlekke golf met een uitgaande bolgolf: 


‚> 
# Edie Lef ì 


| 


wasrin È een konstante, De Fouriergetransformeerde van A volet uit 
ken) 40) 





nà 


Web + í sel lS IEN a 


Í 
î + 
Als we de integratie kunnen uitvoeren is A, dus A, dus ook ï bekend, 


ee " wt EE EN En B 
Hchter, de integrand bevat Y nog, zodat dit onmogelijk is. We bena- 


wt : 
deren Y nu door de uitgaande bolgolf te verwaarlozen; 


zi 


Kek EE en, Tt rid Es 
Alk) = Y Û gE ar) eet en 


Dit is de eerste Bornbenadering. Hogere benaderingen krijgen we, door 
het resultaat uit de eerste benadering te substitueren, etc. Het 
blijkt (zic Merzbacher pag. 227) dat deze benadering goed is als de 
energie van de golf groot is ten opzichte van de inhomogene term «. 


Je WKB-benaiering 


Beschouw de solfvergelijking 


Hai) 0 


u 


En 
waarin Kk e= k(x) een ingewikkelde funktie is. Dan is een vlakke golf 
geen oplossing voor Y , Als k echter langzaam varieert over een golf- 


lengte van Y, kunnen we eens proberen voor welke u(x) de oplossing 


vG) eN iele) 


d 


voldoet (men noemt u wel de eikonaal). Na substitutie vinden we 
‚rè Rr 2 Í = VE 

iv u -— (Vu) + | k(x) j = O 

Als Y een vlakke golf is, hebben we uiteraard 


Hu … |G) } 


Dit doet vermoeden dat Te klein zal zijn voor die k, welke weinig 


veranderen over een golflengte van Y, Hiermee vinden we de benadering 


u 
o 


HI 


+ Ì K(Et) at + konst, 


Xx 
Dit is de nulde WEB (Wentzel, Kramers, Brillouin-) benadering. Schrij=- 
ven we de vergelijking voor u als 


> 2 — 2, 
(Vu) = rt + 1 


Be, 


dan volgt de n+1,_de benadering door substitutie van de n-de in het 


rechterlid: 


% 
Yu zt eld + Shan ) 
n+1 == n 
en dus 
A 2 Ne Es 
Ee Ei : Et ; et et a 

UC) zt J LE Get) + 17 u!) f_ dx! + konst. 

Xx 


>= — 
Als u, () ze u, Ge) verwachten we konvergentie. 


maes erste LD LON LE 


De operator Ì (dew.z. 2, Er En) ) gedraagt zich bij Fourier- 


transformatie als volgt: 


Ï { p(k) oikex } = otk (Tp) + dins = 
= BB) elk IE. L geth | 
VA { pk) ox Ì = Vetter AP + otk VAD) = 
Rek ig ze Walreke} 


We zien, dat in boeide gevallen Fouriertransformatie hetzelfde is als 
kad * 
vervansen van V door ik, Analoos kunnen we bij tijd- in plaats van 


à f o) f ' 
ruimtcetransformatie == vervangen door ij . 








ot 
2e WDedatriees 
Noemen we 2 = Gr, Dr Ps E) dan kunnen de MHD-vergelijkingen ge- 
schreven worden als 
35 == FS 7 wi Feet in EFS ps 98 91 is hmmm 
ot = AGDE = Ko) Be LE) És 1) 
waarin 
H H | 
vx 0 B O& à 0 Ee Om 
| p gp np 
ze 
0 Wes 0 ®) O 6) “Eno o 
He 
| O 0) Me ®) Ö 0 9 “hen 
| p 6) 0) Vx O 0 O 0 
-K(g) = s2 8 8 n on 5 n | 
o o o o 6) le) o 0 | 
Ei 
ï 
0 
Hy, mn ©) 9) ©) “Vy Ve ©) | 
! | 
He 0 -H_ 9 o Vv, 0 B Od 
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EH 
Ö zl { 
| LA 0) 0 0) o Te le) à | 
| h Ex. Z 
E hd. 5 e p ling 8 Ep 
| ad 
| o p Ov. 0 0 0 0 
e d | 
Ig) = | o sp o o Di 0 O 6) | 
Ì =H, He 0) O 6) Vy mi ©) | 
| 
fe) e) e) 0 O 6) 0 o | 
| fe) EH, HE, 0 0 6) =V‚ vy 
H‚ E 
Í V‚ 0 0) O O ne o O 
EL 6 | 
| iz Í 
| 9) Vv, 9) o 0 Ei an 9 | 
1 En | 
e) ©) B U D hp lins 0 | 
0 o pv, 0 0 0 0 | 
-M(g) = | o le) 5 0 Vv o O o | 
| -H‚ o H, © 0 Vz Ö “Vs | 
| | 
| 0 En, EH, 0 0 0 Vs "Vv, 
! ï 
| 6) 0) 0) 0) 0) e) o o / 
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Omdat deze metrices geen differentiaties meer bevatten, proberen 


we ze te diagonaliseren. De karakteristieke vergelijkinsen zijn res- 


pektievelijk 

4 ke 

mjm en se ade BS PL eg 
be lt a ax ax a ax J 
6 L kb 

dh a (s +8” + C Sh + no Ea £ + e°C re Ö 
a S a ay ay a ay J 
"6 Li @ 2 è è 2 Ë è bie} 

À _ G H De) 
aA, La q (s” + Bt De ) + q bed ts keld O0 


waarin Ans Ans Az de eigenwaarden van K, L en M, en ma = Pt A, n 
q 


nk An f A Az . Zoals gewoonlijk is 


[aa 


> 


on en 
\ 


We vinden direkt zes eigenwaarden: 





Dn Es Ser ee ) 


Ú 
Î 


a 
Ar = An = Asy = 0 


Az is TV Aap De Ty Azo pe "V, 





eh, 


Om de overige eigenwaarden te vinden, moeten we een derdegraads=- 
ed : è è è ; ERGE 
vergelijking in m,n resp. q oplossen. Iedere vergelijking van de 


vorm 


U SR + V Us = EU + e Vv Uz = EH + EV 
waarin 
2 4 
e= ej? ali Ü Edw ep? 8 1 
ne Toke 2 ) 
Eb) 7 
% „è 
1 2 1 
n= \ - 5 + ( EW +55 p ) } 
e ĳ 
a ee 
a AE TS 5 


(zie Behnke e‚a. I, pag. 83). In principe zijn de matrices dus te 
diagonaliseren, maar ze hebben niet dezelfde diagonaalvorm, Wel hebben 
alle eigenwaarden het karakter van een snelheid, zodat we in de alge= 


mene formulering K‚ L en M als snelheden mogen interpreteren, 


Dikwijls is een komplexe funktie in een punt niet gedefiniëerd 
omdat de funktiewaarde daar oneindig is. Indien men om zo'n punt heen | 
een omgeving kan vinden, die ‘geheel in het definitiegebied van de 
funktie ligt (met uitzondering van het bewuste punt), dan noemt men 
dat punt een geïsoleerde singulariteit of kortweg singulariteit, Als 
men de funktie moet integreren over een kromme, die singulariteiten 
bevat, komt men in moeilijkheden, Voor een bepaalde klasse singulari= 
teiten, de zgn. polen, is integratie nog mogelijk. 

Gegeven een singulariteit a. Dan kan een komplexe funktie f in 
een ringvormig gebied om a worden ontwikkeld in de Laurentreeks: 


+0 d 
fla) == & a, (z-a) 
—oo 
mits het ringgebied geen singulariteiten bevat. Voor de koëfficiënten 


onderscheiden we drie gevallen: 


1°, Voor alle n <0 is Nn O ‚ Dan noemt men a een ophefbare sin- 


gulariteit; men kan nl. definiëren f(a) = ag 


2". Slechts eindig vele a, met n < 0 zijn ongelijk nul; er is dus 
een geheel getal N waarvoor an Zo en a, = 0 voor n<N, 
Dan heet a een pool van f met de orde N. Het asymptotische ver 


loop van de funktie dicht bij a wordt aangegeven met 





25, 


fla) = O0 ( (aj ) 


of met 
_N 
fla) my (z-a) 


B Oneindig vele a, met n <0 zijn ongelijk nul; dan noemt men a 


een essentiële singulariteit. 


We beschouwen nu funkties die overal analytisch zijn, dew.Ze 
ontwikkelbaar in een machtreeks, uitgezonderd eventueel ter plaatse 
van polen, Deze noemt men meromorfe-funkties. De koëfficiënten a, van 
de Laurentreeks van zo'n funktie berekent men met de integraalformule 
van Cauchy (Zie Morse and Feshbach, pag. 355): 


dj 
@ni 


í =n- 
a, = | £G) Goa) 71 
Y 
waarin y een gesloten gladde kromme om a, die geen singulariteiten be- 
vat, Deze regel stelt ons in staat, de kringintegraal van een funktie 
om een singulariteit heen te bepalen. We zien namelijk, dat voor 


n= 1 


i 


' 
as zut j €) 43 
Y 


Deze integraal is dus bepaald door één term van de Laurentreeks, het 


zgan. residu ag Voor iedere meromorfe funktie f geldt de residuen- 


stelling 


} r 
í 
| f(z) dz = 2ni X Res (Ep) 


5 1 
Y 


waarin / een willekeurige gladde gesloten kromme in het definitiege- 
bied van f£, P; de j-de pool binnen y en Res (E‚p‚) het residu van f 
bij de j-de pool. Het aantal polen r hoeft niet eindig te zijn, Deze 
stelling bewijst men eenvoudig uit de voorgaande vergelijking door vy 
op te splitsen in r delen die elk één singulariteit bevatten. 

Met behulp van het bovenstaande kunnen we funkties over polen 
heen integreren, Laat gegeven zijn een reële funktie 


f(x) = a) 


Xa 


waarin g analytisch is. Deze f heeft een pool in a. Stel dat we willen 


weten 
+00 
f(x) 
J A dx 
- 


dan gaan we als volgt ‘te werk. Konstrueer een funktie f(z), waarin z 





wed 
te 


een kouplexe variabele, die voor reële z gelijk is aan f(x). Beschouw 
dan de volgende kromme: 


Hoem deze kromme Y … Dan is 


Zevlak el 
| f(z) dz e= Ard Res (Ea) 
s 
A De bijdrage aan de integraal van het 
% 
ge kleine cirkeltje verdwijnt als we de 





straal voldoende klein nemen; als f zo= 
Ë. 


danig is, dat ĳ f(z) dz naar nul gaat 
over de grote cirkel als we de straal voldoende groot nemen, is 
+00 
| f(x) dk = 2ni Res (f‚a) 
Als het residu bekend is, is de integraal gevonden, In de buurt van a 
gedraagt f zich als 


fla) wo B 4 (z-a)"| 


Za 


dus het residu is B gla) . Dit kan men ook exakt bewijzen. 
(Horse and Feshbach è 4,5), 


De singulariteiten die bij Green-funkties het meest optreden heb- 
ben de vorm 


klad es 
ne _ a 


Er sijn dus twee eerste=orde polen, in x= a resp. X = =a 


‚ De uit- 
komst van de integraal 


oo 
r 
‘ h(x) dx = ) hz) dz 


hangt van de vorm van v af; waar de integraal bestaat hangt af van het 
asymptotische gedrag van g(z) voor z — + 0. 


Er zijn drie niet-triviale krommen Yv: 


|. 
Z-vlak Z-vlak 
Pd 
f 
4 û \ 
! „a 
De hierbij behorende residuen zijn res- 


am 
a 
à en 1 
Ir e bd Es men 
an pektievelijk: (a) . 5E 
pa 





/ | la) 
a 2 de Áo mtedee lak | 
X 


Er zijn voor deze h(x) dus drie typen in- 


: — + o 
‚ tegralen, die we resp. G , G en G° noemen. 
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